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A long-lived coherent state and non-linear interaction have been experimentally demonstrated for 
the vibrational mode of a trapped ion. We propose an implementation of quantum computation using 
coherent states of the vibrational modes of trapped ions. Differently from earlier experiments, we 
consider a far-off resonance for the interaction between external fields and the ion in a bidimensional 
trap. By appropriate choices of the detunings between the external fields, the adiabatic elimination 
of the ionic excited level from the Hamiltonian of the system allows for beam splitting between 
orthogonal vibrational modes, production of coherent states and non-linear interactions of various 
kinds. In particular, this model enables the generation of the four coherent Bell states. Furthermore, 
all the necessary operations for quantum computation such as preparation of qubits, one-qubit and 
controlled two-qubit operations, are possible. The detection of the state of a vibrational mode in 
a Bell state is made possible by the combination of resonant and off-resonant interactions between 
the ion and some external fields. We show that our read-out scheme provides highly efficient 
discrimination between all the four Bell states. We extend this to a quantum register composed of 
many individually trapped ions. In this case, operations on two remote qubits are possible through a 
cavity mode. We emphasize that our remote-qubit operation scheme does not require a high quality 
factor resonator: the cavity field acts as a catalyst for the gate operation. 



I. INTRODUCTION 

Outstanding theoretical and experimental advances 
have been reported in the field of photonic quantum in- 
formation processing, ranging from the experimental re- 
alization of the quantum teleportation protocol pj to 
proposals for quantum error correction and quan- 
tum computation |3(. It has been shown that universal 
continuous-variable quantum computation can be per- 
formed using linear optics (including squeezing), homo- 
dyne detection and non-linearities realized by photon- 
counting positive- valued- measurement Recently, a 
method to implement efficient universal computation 
based on coherent states of light has been suggested and 
shown to be robust against detection inefficiencies |{|. 

As pointed out in the Los Alamos Roadmap for quan- 
tum computing 6], using coherent states of a boson as 
logical qubits is one of the promising ways to realize quan- 
tum computation. However, one of the practical difficul- 
ties encountered in a scheme for coherent quantum com- 
putation is the requirement of a strong Kerr non-linear 
interaction to produce a superposition of coherent states. 
Currently available non-linear dielectrics, unfortunately, 
offer too low rates of non-linearity with exceedingly high 
absorption of the incoming field. In this context, some 
recent proposals for giant Kerr non-linear interaction ex- 
ploiting electromagnetic induced transparency remains 
to be proved to work in the quantum domain 

In this paper, we propose to implement coherent quan- 
tum computation using vibrational modes of trapped 
ions. Since the early days of the quantum manip- 
ulation of vibrational modes for trapped ions, it has 
been clear that strong non-linear evolutions can be ef- 



ficiently engineered using two or three-level ions (in a 
A configuration) interacting with properly tuned laser 
pulses 0, E3, El E3- Furthermore, a long-lived coher- 
ent state has been experimentally reported 0, . This 
opens a way to the exploitation of vibrational states as 
the elements of a quantum register in a quantum pro- 
cessor. For the purposes of scalability, arbitrarily large 
quantum registers have to be considered. One way is 
to work with a chain of ions in the same trap, exploit- 
ing not just the vibrational modes of the centre of mass 
(CM) but the collective vibrational excitations of the 
chain (see 0,0] and references within). Another way to 
realize the scalability is to take advantage of the recently 
demonstrated coupling between cavities and single-ion 
traps 0, 0| , which is the scheme used in this pro- 
posal. In our architecture, an array of many individually 
trapped ions constitute the quantum register. The local 
processors are interconnected via an effective all-optical 
bus realized by a cavity mode coupled to the different 
traps. We will not require a perfect cavity for our pro- 
tocol and the cavity field mode never becomes entangled 
with the ions of the register (the coupling between two 
different ions being realized via a second-order interac- 
tion only virtually mediated by the cavity field). 

In this paper, we also propose an efficient discrimina- 
tion of the four quasi-Bell states embodied by entangled 
coherent states In this respect, our detection scheme 
does not require the complete map of a quasi-Bell state 
onto the discrete electronic Hilbert space of the trapped 
ions |l5j| . The detection is performed exploiting the ad- 
ditional degree of freedom of the vibrational states rep- 
resented by their even- and odd-number parities. It is 
worth stressing that the Bell-state discrimination can be 
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accomplished, in our set-up, both locally (exploiting two 
orthogonal vibrational modes of a single trapped ion) and 
remotely, where the Bell state is the joint state of two vi- 
brational modes of a linear two-ion crystal. 

The paper is organized as follows. In Section [H] we 
describe the coupling scheme used in our proposal and 
address the issue of the preparation and single-qubit ma- 
nipulation of coherent states. In this context, the gen- 
eration of even/odd coherent states and entangled co- 
herent state is discussed. We perform some quantitative 
investigations to prove that this coupling scheme allows 
for highly efficient quantum state engineering. In Sec- 
tion II I II we propose the architecture for a distributed 
quantum register of many individually trapped ions in- 
terconnected by a cavity held mode. This proposal al- 
lows for vibrational quantum state transfer between two 
remote ions. We quantitatively address a non-trivial ex- 
ample. Section IIVI is devoted to the description of a 
scheme for almost complete Bell-state measurements per- 
formed combining vibrational-mode manipulations and 
electronic-state detection. The ability to achieve a high- 
efficiency discrimination of the four coherent Bell states 
is exploited. In Section we describe how to realize an 
entangling two-qubit gate that, together with the single- 
qubit rotations, allow for universal coherent quantum 
computation. 



II. HAMILTONIAN FOR QUANTUM STATE 
ENGINEERING 



The system we consider is a two-level ion coupled to 
a bichromatic held, detuned from the atomic transition. 
The trap tightly confines the ion in the x — y plane (as 
sketched in Fig. H ( a ))- The energy scheme is shown in 
FigH(b). The external fields (treated classically) illu- 
minate the ion in opposite directions in the z = plane 
and both can have a component along the x and y axes. 
We assume the trap to be anisotropic, with io x > to y and 
Ax cm (Ay cm ) the ground-state width, in the trapping 
potential, along the x (y) direction. The Hamiltonian of 
the system reads (h = 1 is taken throughout this paper) 



H = 



E 



9i°+ 



(1) 



Here, bj (bj) (j — x, y) are the creation (annihilation) op- 
erators describing the quantized position of the CM of the 
ion in the trap, §i = g i e~ l ^ i ' T ~ lu ' it ~ l ^ i take account of the 
couplings between the ion and the i-th laser (i = 1, 2) of 
its frequency Ui, wave- vector = ki V , 0) and phase 
4>i. Here, f is the vectorial operator of the CM position 
and <7+ = &_ = \e)(g\. The ion's transition frequency is 
labelled by oj eg . In a rotating frame and in the limit of 
large detuning A i > g^ 2 , £i2,7rad> where Ai =u) eg — u)i, 
S12 = cl>2 — iO\ and j ra d the spontaneous decay-rate of 
the ion from |e), the atomic excited state can be adi- 
abatically eliminated. After some lengthy calculations 
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FIG. 1: (a): The physical set-up. A two-level ion is trapped 
in a bidimensional trap (ui y < u x ). Two monochromatic laser 
field (classical) excite the ion. (b): The energy level scheme. 
The trap is quasi-harmonic so that, in the limit of resolved 
sidebands, the different excitations of the x and y vibrational 
degrees of freedom are spaced as shown. The lasers are off- 
resonance. 



and using the Campbell-Baker-Haussdorff theorem, the 
Hamiltonian can be written as 



H 



2^ u i b i h i 



r, r, °° v /"„•„> \n+m 

Ai 2— 1 2— 1 n\m\ 

n,rn—0 p.q—0 



(2) 



p\q\ 



-iSi2t—i4 



h.c. 



where <f> = cf>\ — <p2 and a numerical factor arising from 
the normal ordering has been absorbed in the Rabi fre- 
quencies gi,2- Here, r]' x = Ax cm Ak x and rj' y = Ay cm Ak y 
arc the effective Lamb-Dicke parameters for the x (y) 
motion respectively and Ak X}V are the projections of 
ki — k 2 in the z = plane. We have neglected the laser- 
intensity dependent ac-Stark shifts due to the disper- 
sive couplings. These energy terms in the Hamiltonian 
can be controlled by stabilizing the laser beams and for- 
mally eliminated by redefining the ground state energy. 
A scheme to cancel the ac-Stark shifts u sing an addi- 
tional laser has been demonstrated in ref. [16j. Properly 
directing the laser beams we can arrange a coupling be- 
tween the two vibrational modes as well as enginee ring 
a single-mode Hamiltonian (when /i or v is zero ) 111. 
In this latter case, if not explicitly specified, we will al- 
ways consider the states of the x mode to embody the 
qubits, while the y mode will be used as an ancilla. 
An interesting feature of the model in Eq. J3J) is the 
possibility to select stationary terms from the Hamilto- 
nian simply by tuning the laser fields to an appropriate 
sideband resonance of the trapped ion's spectrum. In- 
deed, in the interaction picture, the term depending on 
exp [i{s x uj x + SyLOy — <5i2)i] (and its hermitian conjugate) 
appears in H, where s x — (n — m), s y — (p — q). Tun- 
ing Ji2, which excites the proper sideband of the energy- 
level scheme shown in Fig.n(b), we single out stationary 
terms in Eq. J5J), which we want to be dominant over the 
contributions of the other oscillating terms. 

A remarkable range of evolutions is covered by this 
coupling scheme and some of them are particularly rele- 
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vant for the purpose of coherent quantum computation. 
We note that, in the protocol proposed in y|, the lead- 
ing ingredients are represented by the ability to gener- 
ate coherent states and their macroscopic superpositions 
(Schrddinger cat states) as well as multi-mode entangled 
coherent states. To manipulate the states of the ele- 
ments of a quantum register, on the other hand, ref. p| 
prescribes the use of reliable beam splitter operations, 
phase shifts and displacement operations (these latter ef- 
fectively perform rotations in the computational basis). 
A beam splitter (BS) operation has been described in 
ref. an d we need to give details about the engineer- 
ing of the other operations with our model. 

We need a reliable way to generate a coherent state 
of motion in order to work in a computational space 
spanned by the coherent states {|a),|— a)} (which are 
quasi-orthogonal for sufficiently large a € C). Many 
different ways to achieve this have been suggested 
Here, we note that, if 612 = lu x (i-e. if we excite the 
first red sideband of the x motion) and the two fields 
have no projection onto the y axis, the stationary term 
H cs ~ (igig2Vx/ Ai) b' x + h.c is selected, assuming <fi = 0. 
This energy term gives rise to a unitary evolution that 
corresponds to a displacement D x (a) = exp(ab x — a*b x ) 
in phase-space |5j, ll7j a — g\g^r]' x tj Ai and t the in- 
teraction time. However, if we want to give an esti- 
mate of the accuracy of this state engineering procedure, 
the effect of the non-stationary terms in the Hamilto- 
nian has to be quantitatively addressed. In order to 
do it, we consider the formal relationship between our 
coupling scheme and the system in ref. where the 
three-level A configuration can be mapped onto our own 
coupling-scheme when the adiabatic elimination of the 
excited state of the ion is performed. It is shown in 
that by considering an anisotropic trap with a suffi- 
ciently large ratio lo x /lo v , allows us to neglect additional 
accidental resonances in Eq. © . To generate a coherent 
state, we estimate lo x > 3uj y is enough. On the other 
hand, the coupling factors relative to the non-stationary 
terms in Eq. @ are sensibly smaller than the rate at 
which the coherent state is generated once we guarantee 
id/VxAi) <C 7 with 7 = to x /g a dimensionless parameter 
which, experimentally, can be 7 > 20. For the sake of 
definiteness, we have assumed gi — 92 = .9 and, given 
that S12 <C ^1.2, we have taken |fci| ~ |&2| = k (so that 
i]' x ~ 2rj x ). For the realistic value r\ x — kAx cm — 0.4 , 
the above constraints require Ai 3> <?/5. 

For our quantitative analysis, we choose lu x = 4L0 y and 
Ai = hg. Then, we retain the terms, in H, which oscillate 
at the frequencies w x , uj v , w x ±ujy so that the dynamic 
generator we consider is H trU e = H cs + H nst , where H nst 
collects the non-stationary terms we discussed. With the 
above choices for the relevant parameters, we look for the 
overlap |(a|^) t | between the coherent state \a) we want 
to generate and the state \tp) t = exp {— iH true t\ \$) x . Nu- 
merically, we are limited by the dimension of the com- 
putational space. We thus take \a\ = 1 and truncate the 
basis to {|0) , .., |5)} x , with \n) x (n — 0, ..,5) indicating 
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FIG. 2: (a): The overlap between a coherent state \a = 1) 
and the state generated by the non-ideal Hamiltonian Htrue 
that contains non-stationary terms at low frequencies up to 
second order in the Lamb-Dicke parameter. In this simula- 
tion, r\ x = 0.4. Other parameters are as in the text, (b): 
Comparison between the perfect overlap (solid line) and the 
one obtained with r\ x = 1.0 and 7 = 10 (dashed line). 



phonon-number states. On the other hand, the state of 
the y mode should not be affected at all by the desired 
evolution. Assuming |0) for the initial state, it is reason- 
able to take {|0) , |1)} for the evolution of this state due 
to Htrue ■ The motional state of the ion at time t is there- 
fore expanded as \tp) t = J2n,oJ2 m ,o A 

nm[t) \n, m) X y, 

with E n ,oEm,o \ A n m \ 2 = F The overlap reads 0(t) = 
■^ELo^noW/^' (with Af a normalization factor) 
that can be evaluated once we solve the set of differential 
equations obtained projecting the Schrddinger equation 
for \ip) t onto the |n, m) states. The result is shown 
in Fig. [21 (a). The overlap becomes perfect when the 
rescaled interaction time is gt = Ai/ (gn x ) = 6.25. Fur- 
thermore, we have checked that with the above values the 
efficiency of the process is insensitive to an increase of the 
Lamb-Dicke parameter up to r] x ~ 0.9. For a larger rj Xl 
some deviation from the ideal case is observed. Another 
parameter that is relevant in this investigation is 7. Re- 
ducing it means lowering the oscillation frequencies of the 
non-stationary terms in H true ■ This spoils the efficiency 
of the entire state-engineering process. An example of 
this effect is given in Fig. [3 (b). It is worth stressing 
that, even though the generation of a coherent state with 
just a small amplitude a has been checked here, this will 
also apply for a larger (in principle arbitrary) amplitude. 

For a single qubit operation, we first consider the ro- 
tation U z (8/2) around the z axis of the Bloch sphere for 
the qubit {|±a)} . This rotation is very well approxi- 
mated by displacement operation D x (ie), with e = 6 /4a 

(0e [0,2*]) i- 



UZ[ 2 



e l 2 
e~ ik 



(3) 



e 2iea e (a+ie)bl —(a—ie)b x 



Indeed, D(ie) \a) x 
e 2lae \a + ie) x — e 1 ? \a) x , where the condition a ^> e has 
been assumed (keeping the product at always finite). In 
practice, for a > 2, a small e is sufficient to perform a 2-7T 
rotation. 
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The rotation U x (±w/4) @ around the x axis of the Overlap 
qubit's Bloch sphere can be performed by a Kerr in- 
teraction H K = XK(bib x ) 2 . We now briefly describe 
the procedure to obtain such a Hamiltonian. Align- 
ing the laser beams along the x axis (to get /i = 
1, v = 0) and arranging their relative detuning 5%2 = 0, 
we select a stationary term proportional to b x b x in 
Eq. J2J), with the corresponding rate of non-linearity 
Xk = 2gi52?7z 4 /Ai- We exploit the canonical commuta- 
tion rules between b x and b x and the relation [n x , n x ] = 
(n x = b x b x ), to obtain exp {-iH K t} \ x t=n/2 \Tia) x = 
(\a) ± i |— a)) x @- This macroscopic superposition of 
coherent states is the result of the rotation in the Bloch 
sphere {| ± a)}. 

In order to check the effects of the non-resonant terms 
in the Hamiltonian of our Kerr non-linear evolution, we 
have conducted an analysis similar to the one previ- 
ously performed to generate a coherent state. This time, 
H t rue = Hk + H nst contains terms oscillating at the 
lowest frequencies and up to the fourth power in r\ x . 
The stationary term Hk dominates over H nst because 
the effect of the high-frequency oscillating terms is av- 
eraged out from the effective dynamical evolution of the 
qubit. We retain the same values used before for the 
relevant parameters in our calculations, showing that 
they are suitable for this effective rotation too. Our 
model is sufficiently flexible not to require further ad- 
justments of the set-up. We consider the transformation 
\a = -i) x -> \catl) x = -^{\a = l)+i\a=-l)) x and 

use the truncated phonon-number basis {|0) , .., |5)}® 2 to 
solve the projected Schrodinger equations 



ixy(p,q\d t \i>(t)) xy = xy (p,q\Htrue\lf){t)} xy , (4) 

with \ip(0)) xy = \—i) x |0) y and the decomposition 



m)) xy 



A 



(t) \nm) 



5 



The normalization 

2 



of the wave-function implies m=0 |-A nm (i)|' = 1 
as before. In Fig. [3] we compare the ideal overlap 



O id (t) = \ x (catl 
overlap O trU e(t) = 



-iH K t 



| — i) | (dashed line) to the 
(solid 



1/ <catl > 0|e- iff *~« t |-i > 0) 9il) 

curve), whose time behavior depends on the coefficients 
Anm{t), for counter-propagating lasers with g± t 2 = g 
and r] x — 0.4. It is apparent that the dynamical evolu- 
tion governed by H true leads to the desired superposition 
state at the expected rescaled time gt = 5ir/2r]' x ~ 38. 
The match between the two curves is very good up to 
the rescaled interaction times we show. Being a rota- 
tion in the space spanned by the coherent states with 
amplitude a = ±1, the behavior of the overlap is pe- 
riodic and replicates itself for interaction times larger 
than those shown in Fig. [3] This numerical simulation 
suggests that Kerr non-linear interactions and, thus, the 
U x (±n/4:) rotation, can actually be performed quite ef- 
ficiently in this set-up. We estimate \k — 20tt KHz, for 
g = 5007T KHz 0, ^| , corresponding to an interaction 
time t = it/2xk — 25 /isec. On the other hand, the effec- 
tive life-time of the excited state |e), in the non- resonant 




FIG. 3: Comparison between the perfect overlap Oid(t) (bro- 
ken line) and the true one Otrue(t) obtained numerically solv- 
ing the Schrodinger equation governed by Htrue- In these cal- 
culations, g x = 0.4 and g = 500-7T KHz. Other parameters are 
as described in the body of the paper. 



regime we are considering, is T spont = Al/(g 2 j rad ) ~ 
1 sec if we use metastable levels of an optical transi- 
tion (for example the S^n — > -D5/2 transition in 40 Cet + 
whose excited level has a natural life-time of about 1 sec 
( lrad ~ lsec) dill El). 

We need here to make a remark: the analysis we have 
performed always assumes that the ancillary y mode is 
prepared in the vacuum state. This is just for mathemat- 
ical convenience. We have derived the equations of mo- 
tion for the case of an initial coherent state \(3 = 1) of the 
y-motional mode. Here, again, a small amplitude of the 
coherent state is taken because it is then possible to trun- 
cate the computational phonon-number basis, consider- 
ably simplifying the calculations. We have concluded 
that the comparison between Ovid(t) and Ovt rU e{t) shows 
the same qualitative features seen in Fig. [3] We conjec- 
ture that the same conclusion holds regardless of the state 
in which the y degree of motion has been prepared, if the 
parameter in H trU e are kept within the range of validity 
of the approximations above. 

With arbitrary rotations around the z axis (imple- 
mented via effective displacements) and 7r/4 rotations 
around the x axis, it is actually possible to build up 
any desired rotation around the y axis of the Bloch 
sphere. This, in turn, allows us to arbitrarily rotate 
the qubit around the x axis ,5]. The two operations 
we have demonstrated are thus sufficient to realize any 
desired one-qubit rotation. In particular, the sequence 
U z (it I A)U X (it I '4)U Z (n j '4) \±a) x realizes the transforma- 
tion |±a) x — > (1/V2) (\ot) x ± |— a) x ) that is a Hadamard 
gate. The resulting states, here, are the the so-called 
even (for + sign) and odd (for — sign) coherent states as 
they are the superposition of just even and odd phonon- 
number states, respectively 0- An interesting feature 
that will be exploited later is that even and odd coher- 
ent states are eigenstates of the parity operator (— 1)" 
(n is the phonon-number operator) with eigenvalue ±1, 
respectively. We will discuss later the role these states 
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have in coherent quantum computation. 

As a final relevant case treated here, we now con- 
sider the engineering of a bimodal non-linear interaction 
suitable for the generation of entangled coherent states 
(ECS) 0. This class of states will be represented as 



\4>±)=N£{\a,a)±\- a) -a)} 



\4,±) = N£{\a,-a)±\- 



(5) 



States \4>±) and \ip±) can be generated by superimposing, 
at a 50 : 50 BS, a zero-phonon state |0) with an even and 
odd coherent state, respectively. As an example, suppose 
that, via the procedure described above, we have created 
an even coherent state of the x motional mode while y 
is in its vacuum state. Arran ging a BS interaction be- 
tween x and y phonon modes 1 1 1| . the joint state of the 
two vibrational modes is then transformed into one of the 
entangled coherent states above. Alternatively, \<f>±) and 
\ip±) can be produced using the cross-phase modulation 

with Xcp the rate of non- 
, |/3) , this interaction pro- 



Hamiltonian H cp = Xcpn x n y 



\a) x + \-0)) v + \~a) x (1/5) - \-0)) y 
0, which can be reduced to the form of 



linearity 0. Starting from \a 
duces \ecs) xy a 
when Xcpt = i" 

ECS in Eq. (J5J through singlc-qubit manipulation. Thus, 
having already discussed how to perform one-qubit op- 
erations, we concentrate here on the generation of this 
kind of generalized ECS. 

By inspection of H cp , we recognize the necessity of 
an interaction symmetry in the two vibrational modes. 
This can be obtained directing the lasers at 45 and 
225 degrees with respect to the x axis. In this case, 
6 12 =0 has to be set in order to select the stationary 
term XcpO x b x b y b y in Eq. © with Xcp = 2g 1 g 2 r]' x 2 'q' y 2 / A 1 . 
The other terms in the Hamiltonian are rapidly oscillat- 
ing and negligible if the same dynamical conditions we 
commented above are assumed. We consider O ecs (t) = 



zy (&CSid I 



-iH tr 



|a = i,/? = i; 



■nj 



with Hi , 



= H, 



cp 



H nst the Hamiltonian containing both the desired non- 
linear interaction H cp and all the relevant non stationary 
terms. We have taken \ecs id ) xy cx ll)^ (|1) + |— 1))^ + 
I — 1)^ (|1) — I — X)) , where all the states appearing in 
this expression are coherent states of their amplitude 
M = 1/3 1 = 1- To evaluate O ecs (t), the computational 
basis has been truncated to {|0) , .., |5)}® 2 , as usual. The 
results are shown in Fig. The dashed line represents 
the ideal behavior of the overlap, that is, its time depen- 
dency when just the ideal interaction H cp is considered. 
This curve is contrasted with the overlap obtained when 
the full Hamiltonian H trU e is taken. The mismatches 
between the curves are very small and the overall com- 
parison is excellent. Here, r\ x = i]y = 0-4 and all the other 
relevant parameters are the same as in the previous sim- 
ulations. The rescaled time gt ~ 77, where the overlap 
is almost perfect, is equivalent to an effective interaction 
time of t ~ 50 [isec, using the same parameters of the 
previous calculations. The scheme appears, thus, to be 
robust against the spoiling effects of the non-stationary 



terms and is efficient within the coherence times of the 
physical system we consider. 

Overlap 




FIG. 4: Comparison between the perfect overlap (dashed line) 
and the true one £ca {i) obtained solving the Schrodinger 



equation governed by Hp, 



In this simulation, rj x 



0.4 



and 171,2 = g = 5007T KHz. Other parameters as described in 
the body of the paper. 



III. COUPLING BETWEEN MOTIONAL 
DEGREES OF FREEDOM OF INDIVIDUALLY 
TRAPPED IONS 

So far, our discussions have been limited to the case of 
a single ion. Unfortunately, considering the vibrational 
modes of just a single trapped ion is a substantial lim- 
itation on the computational capabilities of our device. 
However, we can take an advantage of some recent exper- 
iments that demonstrate the coupling between trapped 
ions and a high-finesse optical cavity to design the 

register of a vibrational quantum computer as formed by 
several remote and independently trapped ions. Here, we 
describe in detail a mechanism to couple the motional de- 
grees of freedom of different elements of such a quantum 
register. In the experiments reported in refs. 0,0, a 
coherent interaction is established between an ion and a 
cavity mode. The optical transition between metastable 
states of a Ca + ion has been used to embody a qubit 
and the atomic responses to both temporal and spatial 
variations of the coupling have been analyzed . These 
impressive experimental achievements justify and moti- 
vate a posteriori the two-level model in the optical range 
of frequency we have assumed here. The system sketched 
in Fig. [5J is based on the linear geometry of the ion 
trap-optical cavity interfaces demonstrated in 0, 0|. 
This set-up is suitable to store a linear ion crystal repre- 
sented by a row of aligned traps, which are mutually in- 
dependent and spatially well-separated. The cavity field 
mode (here assumed to be a transverse TEMoo mode of 
a near-confocal resonator) is described by its bosonic an- 
nihilation (creation) operator a (aJ) and is aligned with 
the x-axis of the bidimensional traps. The interaction 
with each ionic transition is off-resonant with detuning 
Aj (i = 1,2) respectively (assumed to be different for 
sake of generality. The mathematical approach is simpli- 
fied if Ai = A2). Two external fields, Eli, -E-L2 excite 
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1 and 2), reads 



FIG. 5: The set-up for a motional state transfer procedure. 
Two trapped ions are placed along the x axis of a mono-mode 
optical cavity having a standing-wave configuration. The ions 
are coupled to the cavity field mode and to two external laser 
pulses. The laser fields illuminate the ions along the y direc- 
tion. This realizes an effective coupling between the y vibra- 
tional modes of the trapped ions. 



the ions and are directed along the j/-axis. As we will 
see, this effectively couples the y modes of the ions. 

We assume a standing-wave configuration for the spa- 
tial distribution of the cavity field with the ions placed 
at the nodes of a cosine function 0, ^| . In a rotating 
frame at the frequency u>li of the laser Eli and in the 
interaction picture with respect to the free energy of the 
resonator, the Hamiltonian of our system reads 



Ai 



j=x,y i=l 



(6) 



Here, S c = uj caV ity - wli and A i2 = ljli - ^L2 while 
£li (i = 1) 2) is the Rabi frequency of the interaction be- 
tween ion i and laser En and <fi is the phase-difference 
between the lasers. The condition Ai^ > £1,1, £1,2, go 
underlies Eqs. 10, where the electronic excited states of 
the ions have been adiabatically eliminated. An intuitive 
picture of the dynamics of the system in this set-up is 
gained in the overdamped-cavity regime (or bad cavity 
limit) , where the cavity decay-rate n is much larger than 
any other rate involved in Eq. JBJ. In this case, the cavity 
mode, which is detuned from the ionic transitions, rep- 
resents an off-resonant bus that is only virtually excited 
by the interactions with the ions and can be eliminated 
from the dynamic of the overall system. In a formal 
way, we can consider the evolution equation of the cavity 
field operator a' and impose that its variations are neg- 
ligible within the time-scale set by the effective coupling 
fij = goSLiVyi/ A-i ( that implies the cavity field mode 
has already reached its stationary state). This results in 
an effective interaction that, in a rotating frame at the 
frequencies of the traps (supposed to be identical for ion 
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i<p 



- 6 T b e 



-i<j> 



(7) 

where the condition A12 = Wj,2 — <*-Vi and the Rotat- 
ing Wave Approximation (RWA) have been used. This 
interaction models a BS operation between motional de- 
grees of freedom belonging to spatially separated trapped 
ions. This interaction is useful for entanglement genera- 
tion and motional state transfer, where the states 
and |0) 2 are swapped, with \ijj) yl being completely ar- 
bitrary. This is exactly what we want to realize for the 
purpose of motional state transfer. If the ion crystal is 
larger than two units, two specific ions can be connected 
by exciting them (and only them) with the laser fields. 
The other trapped ions will be unaffected by the coupling. 
Once the local interaction between the x and y motional 
modes of a specific ion has been performed (according 
to a given quantum computing protocol), then the state 
of the y mode can be properly transferred to another 
ion of the crystal, labelled I, that has been prepared in 
\0) yl . However, for the sake of realism, in what follows 
we pursue the analysis restricted to a two-ion system and 
give some more insight in the process of motional state 
transfer. 

The assumed bad cavity limit is particularly conve- 
nient to isolate the dynamics of the motional modes from 
that of the bus. Indeed, a full picture of the evolution 
of the system is gained by the master equation (in the 
interaction picture) 

d t p = -i [H' ic , p}+ k (2a / oa t - {a f a, p}) = (£ + C cav )p, 

(8) 

with p the total density matrix of the ion 1+ion 2+cavity 
system and, taking 6 C = u> y \, A12 = u) y i — u> y i, it is 

h.c). We have used the no- 
We now go to a dissipative 



tation C p = — i[H' ic ,p). 



pjcture defined by p — e c ^ l p and exploit the relations 
C cav [a, p] = [a,(C cav — k)p\ , C cav (ap) = aC cav p + nap 
(and analogous for aJ) |2(|. After some lengthy calcula- 
tions, Eq. JHJ) reduces to 



dtp = ig^Y.^ 1 {e- K %i,pW + e Kt b yl [a\p] - h.c.} 



(9) 



with C\ (£2) an effective super-operator obtained by col- 
lecting all the terms in Eq. having the e~ Kt (e Kt ) pre- 
factor. To isolate the vibrational degrees of freedom, we 
trace over the cavity mode. We obtain d t p v — £1 (e _Kt p), 
with p v — Tr cav (p). This master equation still in- 
volves the cavity variables because of the presence of p. 
In order to remove these dependencies, we go back to 
Eq. (JSJ), integrate it formally and multiply it by e~ Kt . 
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In the limit of large k, we can invoke the first Born- 
Markov approximation and set p — p v ® p C av,ss with 
Pcav,ss the steady state of the cavity mode. It is, then, 
e- Kt p(t) ~ J °° C 2 (pv(t) <E> Pcav,ss)e~ Kt 'dt' and 



d t pv = Tr cav { Ci ( J C 2 {p v {t) ® Pcav, ss )e K dt' 



E 

2 



Qo^Li^LjVyiVyj 



lyiVyj f 

K I 



i=l 

V^i^2 ^2b yl p v bl 2 + 2b y2 p v b\ x - {fcj^V + 6^1. P«} 

(10) 



with r.; = flf/n (i — 1,2). This is the reduced master 
equation in our study. So far, we have not included the 
relaxation terms due to the decay of the motional ampli- 
tude of the ion modes. However, these can be included in 
the above derivation simply adding to Eq. © the Liou- 
villian terms proportional to the vibrational decay-rate 
7„ (assumed to be equal for the x and y motion). These 
terms do not depend on the cavity operators and are 
left unaffected by the adiabatic elimination of the field 
mode. Their inclusion will eventually result in a modifi- 
cation of just the single mode effective rates according to 
I\ — > Ti + 7„. Of course, there is still much more to un- 
derstand about the mechanisms that lead to vibrational 
decoherence [2l]. However, some estimates of j v put it 
in the range of tens of milliseconds (see @ and Roos et 
al. in |21|) and, as we will see, we estimate them to be 
longer than the effective interaction times required for a 
complete motional state transfer. Thus, from now on, 
we drop 7„ from our analysis. Eq. (|10fl can be projected 
onto the phonon-number basis {\n,m)} yl y2 to give ef- 
fective evolution equations that are used for a numerical 
estimation of the dynamics of p v . 

As an example of motional state transfer, we quanti- 
tatively address the case oi \ip) yl = \/2/5 (|0) — |l)) al + 
y^T/5\2) yl , where (|0), |1), \2)) y i are phonon-number 
states, being prepared in ion 1. Here, the choice has 
been completely arbitrary (any other state could have 
been taken). However, this example offers us the pos- 
sibility to see the influence of our protocol on relative 
phases in general linear superpositions. Furthermore, the 
possible leakage into the Hilbert space complementary to 
the one spanned by our computational basis can be in- 
vestigated. For quantitative calculations, we restrict the 
basis to {|0) , .., |3)}® 2 . For the transfer protocol to be 
effective, the state of the y2 mode must be prepared in 
|0) y2 . Some comments are necessary in order to clarify 
the protocol. From now on, we refer to the ion whose y 
motional state has to be transferred as the transmitter 
while the receiver is the ion prepared in |0) . The inter- 
action channel between the transmitter and the receiver 
is open if and only if both the ions are illuminated by 



the external laser fields. This means that, once one of 
the lasers is turned off, the transfer process stops and 
the interaction channel is interrupted and unable to fur- 
ther affect the joint state of the two ions. On the other 
hand, the effective interaction has to last for a time suf- 
ficient to complete the transfer. We find that tempo- 
rally counter-intuitive laser pulses have to be applied to 
the system of transmitter- receiving ions. In particular, 
an efficient motional state transfer is achieved if the ef- 
fective coupling rate SI2 decreases while fii increases in 
such a way that J T dt' ^T 1 {t')T 2 {t') = = tt/2, where T 
is the total interaction period and time-dependent laser 
pulses have been assumed. An example of such pulses 



-rt 



is given by r x (t) = T 2 (-t) = Te" 14 / [e 

We have assumed Eli = £l 2 = £l, Vyi — r ly2 = Vy an d 
Ai = A 2 = A, so that f = 2g^£l^/A 2 K. The time 
behaviors of Ti, T 2 and y / TiT 2 are shown in Fig. [S] (a). 
To satisfy the conditions for the adiabatic elimination of 
the cavity mode (fii, VL 2 <C k), we have taken T = 0.03. 
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FIG. 6: (a): We plot F^t) (dashed line), F 2 (t) (dot-dashed) 
and s /F 1 (t)F 2 (t) (solid line) as a function of the rescaled in- 
teraction time Kt and for F = 0.03. The interaction channel 
is open at t — — 200k; - . The fidelity does not reach 1 be- 
cause of the effective dissipative part in the evolution of the 
vibrational state that, actually, raises the vacuum component 
of the joint vibrational state, (b): Fidelity of motional state 
transfer |?/>,0) 1 y2 — > |0, ip) yl y2 as a function of the interac- 
tion time nt. M er the interaction channel is closed, the state 
of the two vibrational modes is stationary. 



The time integration of the effective coupling rate 
range where the amplitudes of the laser pulses are non- 
zero gives 1? ~ 7r/2, so we expect that, after the effective 
interaction mediated by the cavity mode, the yl motional 
state has been transferred to y2. This can be seen plot- 
ting the fidelity T^(t) = v \, y2 (0,V>|$(pv) \^^)y\.. y2 -> as a 
function of Kt. Here, $ is the map given by the reduced 
master equation we have derived. It takes the initial den- 
sity matrix p^ = | - 0,O) yl a2 ( - i/',0| into p v {t), the solution 
of Eq. (|10fl . The fidelity turns out to be a function of the 
interaction time and is parametrized by the state we want 
to transfer. The results of our simulation are presented 
in Fig.[S](b). The fidelity is very good, reaching ~ 0.9 for 
t > 200k -1 . We have plotted !F^{t) for interaction times 
larger than T to show that, once the effective coupling is 
turned off, the interaction channel breaks down and the 
state of the ions becomes stationary. It is worth stress- 
ing that the fidelity at the beginning of the interaction is 
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non-zero because of the presence of |00) j y2 in both the 
initial and target states. 

The second important point that has to be addressed 
in order to completely characterize the performance of 
the state transfer protocol is the leakage. We can single 
out two different kinds of leakage. One kind is to states 
such as \i, j) yl y2 (i, j £ {1, 3}) which are the states of the 
computational basis having more than 2 phonons in the 
yl mode and some phonons in the yl mode. The other 
kind of leakage leads to states lying outside the compu- 
tational space. The influence of both these sources of 
error can be contemplated looking at the norm of the fi- 
nal density matrix p v (t) = We have checked that 
Tr(p v (t)) = 1 for all the relevant interaction times, show- 
ing that the influence of highly excited phononic states in 
the vibrational modes can be neglected. This indirectly 
demonstrates that leakage of the latter kind is irrelevant 
and the dynamics of the system is confined in the com- 
putational space we have chosen. On the other hand, by 
considering the quasi-norm J2i=o yi,yi{Qi\ $(pv) \^) y i, y 2^ 
we can check how large the influence of populated yl 
states is in the density matrix p v (i). This is shown in 
Fig. [7| (a) where we can see that, after the transient pe- 

(a) (b) 
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FIG. 7: (a): Fhe quasi-norm J2t Oyl , y2 WHp<p)\0i} yl 
plotted versus the interaction time Kt. This quantity considers 
the contribution to the normalization of the density matrix 
given by the states having just the y2 mode populated, (b): 
The degree of mixedness of Pv(t) versus the interaction time 
Kt. The purity of the state is measured by the linearized 
entropy Sl = (16/15) [l - Tr(pl)] that is zero for a perfectly 
pure state (as the initial one, at Kt — —200) and 1 for a 
statistical mixture. 

riod when the contribution by non-empty states of the 
yl mode is relevant, the steady state of the system is 
perfectly normalized. This means that the final state of 
the two vibrational modes does not contain excitations of 
the transmitter. To complete this analysis, we give some 
insight about the purity of the state we get. From this 
viewpoint, the fidelity is not a good tool because it could 
give the same quantitative results for p v (t) describing a 
pure state or a statistical mixture. On the other hand, 
an easily computable quantity is the linearized entropy 
S L (t) = (16/15) [l-Tr{pl{t))} H. This quantity is 
zero for a pure state and reaches 1 if the state is com- 
pletely mixed. We show a plot of <Sx(t) in Fig. (b). 
The state remains highly pure all along the interaction, 
the small degree of mixedness being due to the dissipa- 
tive nature of the effective evolution of the system arising 



from the bad-cavity regime. 

The fidelity and purity of the state is not perfect be- 
cause of the losses induced on the two-mode vibrational 
subsystem by the dissipative bus. However, the state we 
get with this protocol is nearly optimal. Higher quality 
factors of the optical cavity coupled to the ion traps will 
improve the performances. In this case, indeed, we would 
be able to neglect the dissipative dynamics of the cavity 
mode, making the effective interaction between trans- 
mitter and receiver perfectly unitary. The form of the 
effective Hamiltonian, arising from the coupling scheme, 
will be as in Eq. J7J) with the replacement k — > S c as, 
in this case, the condition S c 3> K has to be used. The 
fidelity of the operation and the purity of the final state 
will be ideal. Unfortunately, this is not the case of the 
experiments performed in |l4|. where the cavity decay- 
rate is k ~ 2tt x 102 KHz and has to be contrasted 
to the ion-cavity coupling go ~ 2ir x 134Hz. These 
cavity parameters, an external laser-ion Rabi frequency 
Sl — 30 KHz with Lamb-Dicke parameter r\ y ~ 0.2 and 
detuning A ~ 50 KHz (larger than a typical laser band- 
width ~ 6KHz) allow us to get an effective coupling 
rate T ~ 0.03. This, in turn, gives interaction times in 
the range of 0.1 msec, smaller than the effective lifetime 
of the ion's excited state. It is noticeable that we do not 
require a high-quality factor cavity to achieve an almost 
perfect motional state transfer, a feature of the protocol 
we have proposed that is certainly important from the 
practical viewpoint. 



IV. QUASI-BELL STATES MEASUREMENT 

We next consider the Bell-state measurement needed 
in the protocol for coherent quantum computation [j| . As 
we will see in Section Bell-measurements can be used 
to construct the teleportation-based CNOT suggested by 
Gottesman and Chuang |24|. In our specific case, the 
quantum channel for the teleportation protocol is em- 
bodied by one of the ECS's in Eq. (J5J. For sufficiently 
large amplitudes of their components, the ECS are quasi- 
orthogonal, carry exactly one ebit of entanglement and 
are usually referred to as quasi-Bell states. A complete 
discrimination of the elements of this class is, thus, fun- 
damental in our scheme. It is worth stressing here the 
well-known no-go theorem demonstrating that a never- 
failing, full Bell-state analyzer can not be realized using 
just linear interactions (e.g. beam splitters and phase 
shifters) [25J. More recently, it has been recognized that 
the introduction of a Kerr non-linear interaction |26| or 
the exploitation of additional degrees of freedom of the 
system employed |27| can be used to fully discriminate 
all four Bell states. However, these schemes are designed 
to work with two-level systems and are not relevant to 
the infinite dimensional case we treat. 

The direct detection of the properties of a vibrational 
state is, in general, a hard task to accomplish. On the 
other hand, detecting the electronic state of an ion (or 



9 



an array of ions) is more straightforward and can be per- 
formed using the quantum jump technique, in which res- 
onance fluorescence from a strongly driven atomic tran- 
sition is detected 0, El EH The presence/absence 
of fluorescence in the driven transition reveals the elec- 
tronic state of the ion. Thus, we need a joint interaction 
that changes the internal degrees of freedom of the ion in 
a way that reflects the state of the vibrational ones. The 
measurement of the electronic state of the ion after the 
interaction, then, will give information on its vibrational 
state. To achieve this goal, we start by considering the 
Hamiltonian obtained by applying a standing-wave laser 
field to the trapped ion. In a rotating frame and with 
A = io eg — lol the detuning between the standing- wave 
and the ion's transition frequency, the interaction reads 

H' = Aa z + ui x P x b x + Qcos[r](bl + S x )](<7_ + <t+), (11) 

where f2 is the corresponding Rabi frequency. In the 
dispersive limit A ^> fi, with A well-away from the res- 
onant vibrational frequency (see Schneider and Milburn 
in [2l]]), we can adiabatically eliminate the excited state 
of the ion and expand cos^&J. + b x )\ in power series, re- 
taining terms up to the second order in rj (Lamb-Dickc 
limit). In the interaction picture and neglecting terms 
oscillating at frequency ±2lo x , we get 



n 



H q nd — 2—r) 2 b x b x a z , 



(12) 



where state-independent energy terms have been omit- 
ted. This Hamiltonian is suitable for quantum non- 
demolition measurements of the motional even/odd co- 
herent states discussed above. Indeed, the evolution op- 
erator U qn d(t) = exp {—iH qn dt) does not change the par- 
ity of an even/odd coherent state but phase-shifts the 
electronic state by an amount depending on the vibra- 
tional state phonon-number. Explicitly: 



Uqnd(t) = coa(x q ndtblb x )l - % sin(xqndtb x b x )a z 
= jx^dtVJ* \g}(g\+ e - lx i" dt ~ b *~ h * \e)(e\ 



(13) 



with Xqnd = 2f2 2 ?7 2 /A. If we set Xqndf — t/2, the elec- 
tronic states will be mutually shifted 7r-out-of-phase. 

Now, let us assume that we prepare the x vibra- 
tional mode of an ion in an even/odd coherent state 
N±(\a) ± I — ck)), while its internal state is \g). Then, we 
apply a 7r/2-pulse tuned on the carrier frequency of the 
ion's spectrum. This particular interaction realizes the 
Hamiltonian H car = Q&+ + h.c. that couples \gn) <-> \en) 
(G being a Rabi frequency). That is, it does not af- 
fect the vibrational state The 7r/2-pulse prepares the 
superposition (l/\/2) (|e) + \g)). The standing-wave de- 
scribed above is then applied and the interaction lasts for 
t = n/2xqnd- This step of the protocol is used to write 
the vibrational state on the internal degrees of freedom 
of the ion. Another carrier- frequency 7r/2-pulse mixes 
up the phase-shifted components of the electronic state 
and, finally, the internal state detection is performed via 



quantum jumps. It is worth stressing that the electronic 
state detection is a true projective measurement (in the 
Von Neumann sense) that is able to tell us if the ion was 
in \g) or not. In this latter case, the vibrational state is 
reconstructed depending on the outcome of this last step. 
The described protocol realizes the transformations 



Af+ (\a) + \-a)) \g) [\ia) + \-ia)) \g) 
tf_(\a)-\-a))\g)^{\ia)-\-ia))\e) 



(14) 



Thus, the different parity of the two vibrational states af- 
fects differently the interference between the components 
of the Fourier-transformed state \g). The discrimination 
between even and odd coherent states can be performed 
with, in principle, high accuracy. Each step in the pro- 
tocol can be, indeed, quite precisely performed if a judi- 
cious choice of the parameters is made. The preparation 
of the electronic state superposition can be done off-line, 
exploiting one of the two laser beams that build up the 
standing-wave and reminding that the effective Hamilto- 
nian in Eq. @ does not affect the electronic variables 
of the ion (the manipulation of the electronic state then 
has no influence on the vibrational states). An estimate 
of the interaction time required to perform U(Tr/2x qn d) 
and to achieve the right phase shift leads to ~ 15 usee 
for Q = 7T x 500 KHz, 77 = 0.2 and A = 10 MHz jig. 

This protocol, which was studied for the cavity quan- 
tum electrodynamic model to detect even and odd par- 
ities of the cavity field 30], is useful for the detection 
scheme for ECS. In particular, let us suppose an ECS 
state of the x and y modes of an ion is subject to a 
50 : 50 BS operation. This will give us one of the output 
modes in an even/odd coherent state, the other being in 
its vacuum. In particular 
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B : 



y = Ni{\V2a)±\-V2a)} x ®\0) y , 
y (J) \^±) xy = N£ |0) x ® {\V2a) ± \-V2a)} , 



(15) 



where B xy (ir/A) is the 50 : 50 BS operator H3- Then, the 
following protocol could be used. We prepare the elec- 
tronic state of the trapped ion in the ground state and ap- 
ply a carrier- frequency 7r/2-pulse to get the electronic su- 
perposition (|e) + \g)) /\/2- Then, we arrange the evolu- 
tion Uq n d{t/2xqnd) for the electronic+x- vibrational sub- 
system. Another carrier-frequency 7r/2-pulse on the ion 
mixes the components of the electronic superpositions. 
A quantum-jump detection reveals the internal state and 
the output is recorded. If the result of the measure- 
ment is I g), the entire protocol is re-applied, this time ar- 
ranging the U qn d(t/2xqnd) evolution of the electronic+y- 
vibrational subsystem. If the result of the first electronic 
detection is instead |e), we use a 7r-carrier pulse that re- 
stores \g), before the protocol is re- applied (this can be 
done with, in principle, 100% of accuracy 0)- The differ- 
ent combinations in which the internal state of the ion is 
found allows us for a partial discrimination between the 
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elements in the ECS class. Denoting \even,iy/2a)j and 
\odd,iy/2a)j (j = x,y), respectively, the even and odd 
coherent states whose components have absolute ampli- 
tude liV^alj one can prove the correspondences shown 
in Tab. D 



Initial state 


I st det. 


2 nd det. 


Final vibrational state 


^"t" ' 9 ) xy ,ion 


Is)- 


Is),- 

t* 3 1 ton 


jewen, i\/2a) ^ <g> |0) v 


\4>-,g) ■ 


|e>. 

I / ion 


\g)- 

l** / ion 


|odd,iV2a) T ® |0} v 


1 J 9 ) xy ,ion 


lff)i„n 


lfl) io n 


10)^ ® \even, i\/2ot) ^ 


IV 1- ' 9 ) xy,ion 


ls)*on 


1 *") ion 


|0) T ® odd, iV2a) v 



TABLE I: We schematically present the protocol for a par- 
tial quasi-Bell states discrimination. The initial state (the 
vibrational x and y modes being prepared in an ECS), the 
outcomes of the 1 st and 2 nd electronic measurement and the 
final vibrational state are displayed. 



It is seen that, while the discrimination between 
\<t>-) xy > \^-) xy and {\<t>+) > \^+)} xy is Perfect, this is not 
the case for the elements of the subset {|</>+) , \ip+)} xy . 
The sequence of detected electronic measurements cor- 
responding to these vibrational states is the same and 
there is no way to distinguish between them, following 
this strategy. However, one can exploit the parity non- 
demolition nature of the above procedure. Even if the 
amplitudes of the components of an even/odd coherent 
state are changed (the amplitude transforming from \a\ 
to |i\Z2a|), the parity eigenvalue of these states is pre- 
served. 

Now, if \4>+) xy is prepared instead of \ip+) xy , we end 
up with mode x being populated while y is in its vibra- 
tional vacuum state. The configuration will be contrary 



if \ip+) xy is prepared. Thus, the key of our procedure 

is the discrimination between 10)^ and |euen, iy/2a) x. 
Let us suppose that, after the application of the previ- 
ous protocol and having found a sequence of two ground 
states as a result of the detection procedure (with the 
prepared vibrational state being totally unknown), we 
apply the displacement operator D{—e) to mode x, e be- 
ing a proper amplitude. If \4>+) xy was the initial state, 
the displacement transforms the state of the x mode into 



i2\/2o 



|n/2a - e) + e- i2y/ * ae \-is/2a - e)\ . Taking 

V J X 

2^/2ae = it/2 and a = 2, then we get e ~ 0.27 -C a. 
With this angle of rotation the even coherent state is 
changed into an approximation of an odd one. Wc 
have flipped the parity of the state. Applying now the 
quasi-Bell state detection protocol, as described above, 
the outcome of the first atomic measurement becomes 
\e) ion . Does it help in distinguishing between \4>+) xy 
and \ip+) xlx2 1 The effect of displacement on the state 
® (|*v2 Q! ) + | — ^^ a )) y (that is the final vibra- 
tional state if \?p+) xy is prepared, as shown in Tab. QJ 
is 1 0)a; ~~ > \~ e ) x - For the amplitudes of the coherent state 



a and the angle of rotation e chosen above, however, it is 
\-e) x ~ 0.96 10)^-0.27 \l) x . Applying again the ECS de- 
tection scheme, just before the first electronic detection, 
we get 



\-e) x <g) even,iV2a) (g> \g) lon 

-» (0.96 10) |0)-iO.27|l) |e)) X)ion ® 



even 



(16) 



The probability to get the electronic output \e) ion is given 
by (0.27) 2 < (0.96) 2 . That is, most of the times (approx- 
imately 92% of the times) we will obtain \g) ion , making 
the discrimination between the two states complete. This 
can be taken as an estimate of the efficiency of the quasi- 
Bell state measurement because, in the present case, the 
efficiency of the detector apparatus (that is the efficiency 
of the quantum jumps technique for electronic state's de- 
tection) can be taken, nominally, 100% Our pro- 
tocol for quasi-Bell state measurements can be adapted 
to the case of two distinct vibrational modes relative to 
remote trapped ions. It could be re-designed, mutatis 
mutandis, using two ions and their x motional modes. 
In this case, the internal degrees of freedom have to be 
detected in parallel and not sequentially and the discrimi- 
nation will be based on the comparison between different 
combinations of them. Considering all the operations 
to be performed before the ion's internal-state detection, 
the overall time required for a complete discrimination of 
ECS should be in the range of hundreds of /zsec, which 
is within the coherence times of the system. 



V. CONTROLLED TWO-QUBIT GATES 

In this Section, we consider controlled two-qubit gates 
to complete our discussion on a possibility of qubit op- 
erations using coherent states of vibrational modes of 
ions. As we have remarked in the previous Section, a 
quasi-Bell state detection (both local and distributed) is 
possible. The teleportation-based scheme for a CNOT, 
indeed, exploits the Bell state measurements to perform 
two different steps. We follow the scheme proposed in [B| 
and consider two three-mode GHZ states, |£)„ a i ai ana - 
I0a3 a4 o5' °^ general bosonic modes a0, .., ah. The joint 
state of modes a0 and a3 is first projected onto the Bell 
basis in order to prepare the (un-normalized) four-mode 
entangled state 



Mane = \ a > a )al,a2 l^+)a4,a5 + I" 



- a )al,a2 IV'+) Q4 ,a5 • 

(17) 



This is then used to realize the CNOT gate as described 
in refs. 0, 0- The GHZ states can be built by using 
beam splitters [ll| and single-mode rotations as those 
already demonstrated above. It has been recently rec- 
ognized (see for example ) that the four-mode state, 
being a complicated step to perform, can be prepared 
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off-line and then used in the protocol for the CNOT only 
when it is needed. Overall, we need 8 vibrational modes 
to implement a single CNOT, six of which are used to 
prepare for \r)) and the remaining two for the control 
and target qubits. In our scheme, however, we use two 
vibrational modes per ion so that we require 4 ions and 
many motional state transfer operations. Even if there is 
no in-principle difficulty in doing this, it is apparent that 
the scheme is experimentally challenging, not just for the 
in situ operations to perform (linear and non-linear cou- 
pling between orthogonal vibrational modes) but for the 
transfer protocol (that is the slow and less efficient part 
of our scheme). However, it is straightforward to extend 
the Hamiltonian model in Eq. J2J) to three orthogonal 
vibrational modes (that is to include the z mode in the 
coupling model) using laser fields having (opposite) pro- 
jection onto the azimuthal axis too. In this way, we will 
be able to exploit three modes per ion, altogether, and 
the steps necessary to create \rj) anc can be performed us- 
ing a two-ion crystal in the optical cavity. The projected 
modes etO and a3, having not been involved in the four- 
mode ancillary state, could be used to embody the target 
and control of a two-qubit gate. A single CNOT gate, 
thus, can be realized with just a two-element register. 

Furthermore, our ability to engineer the Kerr non- 
linearity can be used here to reduce the number of tele- 
porting operations we have to perform. We exploit that 
\ ( t ) +) a i,a 3 and l^+)ai,a? are mutually swapped by the 
cross-parity operator (— l) n «™"j. On the other hand, 
are not affected by this evolution. 

Schematically: 



(18) 



' at, a j 



(- 




)ab = 


\i>+)ab> 


(- 




)ab = 


\<t>+) a b, 


(- 


-l)""™ 1 \<j>_ 


)ab = 


l*->.6> 


(- 


■if*™* |v>_ 


)ab = 


\^-)ab- 



The operator (— 1)"° 



is implemented by the 



cross-phase modulation used to create the ECS from sep- 
arable coherent states, as described in Section |H] Here, 
we are interested in the effect of this unitary operation 
on the class of ECSs. Now, a cross-phase modulation 
between vibrational modes al and a3 is assumed to be 
used so that the state |\/2a) al (|\/2a) + |— v2a)) 3 + 

\-V2a) al (\V2a)-\-V2a)) a3 is generated 0. We 
consider the vibrational modes a2 and a4, each prepared 
in the vacuum state, and realize a 50 : 50 beam split- 
ting in the al + a2 and a3 + a4 subsystems. Finally, the 
(— l)" a3n ° 4 produces a state (not normalized) that, with 
local unitary operations (a single-qubit NOT gate), can 
be written as 



\V 



K a )al,a2 |0+)o3,o4 + I" 



a '- a )al,a2\lP-)a3.a4- 

(19) 



This differs from \r)) anc because it involves 
(correlated to |— a, — a) ol o2 ) instead of \?P+}a3a4- 



l<» + >< 
|c>- 



mode xl 



|t>- 
|0 + >< 



Bs(e) 



a3,a4 

We 



mode x2 

FIG. 8: Scheme for a coherent-state CNOT. Target (\t)) and 
control qubit (|c)) are superimposed at a beam splitter (reflec- 
tivity cos 2 (#/2)). The states of the output modes are, then, 
teleported onto the state of mode xl and x2. The boxes la- 
beled T resume the teleportation operations. 



do not try to reproduce the four-mode entangled chan- 

nel Vflanc - 24 J but W<3 S° 011 Witn W)anc and a PPl y 

the protocol for a teleportation-based two-qubit gate. 
It can be proved by inspection that, in this case, a 
controlled-i&y (Qq- ) gate between the control c and 
the target t qubits is realized (up to single-qubit ro- 
tations to be applied, conditionally on the outcomes 
of the Bell detections). In the computational basis 
{\ot, a) , \a, —a) , |— a, a) , |— a, —a)} ct , this can is repre- 
sented by the block diagonal matrix Q CTh = diag[l, ia y ], 
with 1 the 2x2 identity matrix and a y the y-Pauli ma- 
trix. This gate is non-local and is not locally-equivalent 
to a SWAP gate. It is indeed easy to see that Q ay is 
an entangling gate that transforms the separable state 
(A | a) +B \-a)) c ® (C \a) +D\-a)) t into the entangled 
state A \a) c (C \a)+D \-a)) t +B \-a) c {D \a)-C \-a)) t . 
Thus, together with the single-qubit rotations in Section 
II, this two qubit operation can be used to perform the 
universal quantum computation |33j . Moreover, using 
the criteria of refs. ^S^ 7 this gate turns out to be locally 
equivalent to CNOT. 

On the other hand, an important simplification in 
the realization of a two-qubit gate can be achieved if 
we renounce the four- mode entangled channel \rj) anc (or 
W) anc )- This latter can be replaced by two \<t>+) ECS 
states that are used as quantum channels for the tele- 
portations of the output modes of a BS (reflectivity 
cos 2 (6>/2)) that has superimposed the control and tar- 
get qubits (see Fig.(SJ. For proper choices of 9 = n/Aa 2 , 
after the beam splitter operation and the two teleporta- 
tions, the output modes (xl and x2 in Fig. |SJ are in a 
state that is equivalent to CNOT, up to single-qubit op- 
erations. To achieve this result with a significant proba- 
bility of success, however, the condition 8 2 a 2 <§C 1 has to 
be fulfilled. For example, if we take a = 2, 6 = n/36 has 
to be taken (corresponding to an interaction time with 
the external laser fields of some /isec, for the parameters 
used in this paper) and giving a probability of success 
~ 0.92. Explicitly, the control and target qubits can be 
written on the y modes of the trapped ions 1 and 2. The 
ECSs we need for the teleportations are codified in the 
state of the x and z modes. The beam splitter operation 
between the control and target qubits is, in this case, 
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the sole delocalized operation we need to perform. The 
scheme we have described for motional-state transfer is 
exactly what we need to split the input modes in a dis- 
tributed way. All the other operations do not involve 
coupling between the ions of the register. The Bell-state 
measurement is finally performed involving z and y mo- 
tional modes, projecting the x modes onto a state equiv- 
alent to CNOT \c, t) xl x2 . We conclude our analysis with 
a remark concerning the strategy to follow in order to 
discriminate the logical states of the qubit (namely be- 
tween | a) and |— a)). This can be done locally, following 
ref. using a 50 : 50 BS operation superimposing the 
state of the qubit, codified in the x mode of an ion, to a 
coherent state of the ancillary y mode. Then, with proper 
resonant transitions coupling the internal state of the ion 
to the x and y modes and highly-efficient electronic state 
detection, we can ascertain the state of the qubit 0] . The 
generalization of this procedure to a delocalized situation 
can be done exploiting the results shown in Section ITTT1 

VI. REMARKS 

In this paper we have presented a scheme that, ex- 
ploiting the motional degrees of freedom of individu- 
ally trapped ions, could allow for coherent-state quan- 
tum computation 0, . We have addressed a model 
for quantum engineering based on the use of two non- 
resonant laser pulses. By regulating the direction of the 
lasers along the trap axes and tuning their frequencies 
to excite proper sidebands, we realize various linear and 
non-linear interactions, both for the one-qubit and the 



two-qubit operations. To scale up the dimension of a 
quantum register, we have considered a distributed de- 
sign of the quantum computer, each node of the network 
being a single trapped ion. The interconnections between 
remote nodes are established by a cavity-bus coupled to 
the transition of the selected ion [T^.ll4j. Motional state 
transfer, in this way, is shown to be realizable with good 
fidelity and without the requirement of a high-quality fac- 
tor cavity. Finally, an efficient quasi-Bell state discrimi- 
nation is possible, in this set-up, using unitary rotations 
of the states belonging to the ECS class and inferring the 
parity eigenvalues of superposition of coherent states via 
high-efficiency electronic detections. The accuracy of this 
scheme can be, in principle, arbitrarily near to 100% due 
to the exploitation of the additional degree of freedom 
represented by the electronic state of the ion. This fea- 
ture allows us to circumvent the bottleneck represented 
by the no-go theorem in ref. |25j . We have addressed 
the issues of efficiency and practicality of our proposal 
showing that, singularly taken, each step of the scheme 
is foreseeable with the current state of the art technology, 
the main difficulty, up to date, being represented by the 
sequential combination of them. 
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